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ABSTRACT 

The quest to detect prebiotic molecules in space, notably amino acids, requires an understanding of 
the chemistry involving nitrogen atoms. Hydroxylamine (NH 2 OH) is considered a precursor to the 
amino acid glycine. Although not yet detected, NH 2 OH is considered a likely target of detection 
with ALMA. We report on an experimental investigation of the formation of hydroxylamine on an 
amorphous silicate surface via the oxidation of ammonia. The experimental data are then fed into a 
simulation of the formation of NH 2 OH in dense cloud conditions. On ices at 14 K and with a modest 
activation energy barrier, NH 2 OH is found to be formed with an abundance that never falls below a 
factor 10 with respect to NH 3 . Suggestions of conditions for future observations are provided. 
Keywords: ISM: molecules — ISM: atoms — ISM: abundances — ISM:dust, extinction — Physical 
Data and Processes: astrochemistry 


1. INTRODUCTION 

With molecules detected in space increasing in num¬ 
ber and growing in complexity, one of the main goals in 
astronomy today is to find a link between astrochem¬ 
istry and astrobiology (iWincel et al.l l20doll . Prebiotic 
molecules have been already observed that can be con¬ 
sidered as precursors of complex organic compounds like 
amino acids and prebiotic material. N-bearing molecules 
are necessary precursors of such prebiotic material and 
among them hydroxylamine (NH 2 OH) appears to be one 
of the best ca ndidates. Both exper i mental and theo- 
retical studies (|Blagoievic et al.ll2003t iSnow et ^1200711 
have shown that in its protonated form it can react in 
the gas phase with acetic and propanoic acids to yield 
protonated glycine (the simplest non-chiral amino-acid) 
and protonated alanine (the second simplest and chi¬ 
ral amino-acid). However, the question is still under 
debate whether such a synthesis happens in the gas 
phas e or by surface r eactio ns. Actually, in a recent 
work iBarrientos et al.l (j2012l l dismiss the possibility that 
precursors of glycine could be efficiently produced from 
the reactions of hydroxylamine-derived ions with acetic 
acid in the ga s phase. This is confirmed by recent work 
l|Garrodl[201311 using a new three-phase chemical model 
that fully couples gas-phase, grain-surface and bulk-ice 
chemistry, which shows that glycine can be formed much 
more efficiently on grain surfaces than through gas-phase 
reactions. 

Gurrently, the main problem is the lack of astronom¬ 
ical detection (either in the inter- or circum-stellar me¬ 
dia, ISM and GSM) of both hydroxylamine and glycine 
(NH 2 CH 2 GOOH), and this is obviously correlated to 
their low predicted abundances. Searches for hydrox- 
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ylamine have been done by iPulliam et al.l (1201211 ob¬ 
serving seven bright sources toward IRG-l-10216, Orion 
KL, Orion S, Sgr B2(N), Sgr B2(OH), W3IRS5, and 
W51M using the NRAO 12m telescope on Kitt Peak, 
and by iGharnes^ (j2013l l toward L1157-B1 using the 
CARMA Radio Telescope Array. Glycine of cometary 
origi n has been detected in STARDUST- returned sam¬ 
ples ([Glavin et al.ll200§ : lElsila et al.ll2009f) . but it is un¬ 
known whether it represents pristine interstellar mate¬ 
rial. A chemically related species that is a possible direct 
precursor of glycin e, amino acetonitrile (NH 2 CH 2 GN), 
has been detected (jBelloche et al.lI2OO8II in Sgr B2 (n/ 
Glycine detection in the ISM or GSM is claimed to 
be highly plausible using ALMA to wards bright nearby 
sources with narrow emission lines (|Garrodll2013l l. Such 
sources are luminous high-mass young stellar objects 
and low-mass protostars, also called hot cores and cori- 
nos. They are regions surrounded by dense material (10® 
cm“® for hot cores and 10 ® cm“® for hot corinos) and 
warm (~250 K and 100 K, respectively) close to the core. 
This material is much colder farther away. Deep within 
such dense regions, cosmic-ray (GR) ionized helium, 
He+, breaks up the very stable N 2 molecule, liberating 
N atoms that interact with abundant H atoms to yiel d 
NH 3 , which then sticks to cold grains ([Tielensl 1200511 . 
It is worth noting that the composition of hot-cores re¬ 
sembles that of interstellar ices with high abundances o f 
hydrogenated species as H 2 O and NH 3 (jTielensI 1200511 . 
Hot cores and corinos are in general places where many 
terrestrial-like organic molecules are observed. More re¬ 
cently, radiative transfer calculations of glycine emission 
show that low-mass star-forming regions, in their earli¬ 
est stages of cold pre-stellar cores, may be better region s 
for the detection of glycine (|Jimenez-Serra et al.ll201^ . 
Water detection in the low-mass pre-stellar core L1544 
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l|Caselli et al.ll201^ indicates that a fraction of the grain 
mantles has desorbed into the gas phase. This makes the 
observation of glycine plausible, assuming it co-desorbs 
with water . After the first detectio n of nitrogen hydrides 
in Sgr B2 (jGoicoechea et al.ll2004ll . a complex region en¬ 
compassing a variety of physical conditions and in par¬ 
ticular shocks, the new detection with Herschel/HIFI of 
NH (and ND), NH2, and NH3 towards t he external enve¬ 
lope of the protostar IR-AS 16293-2422 (|Hilv-Blant et al.l 
120101 : iBacmann et al.l[2010f l has not improved our knowl¬ 
edge of the interstella r chemistry of nitrog en hydrides. In 
a recent publication, ILe Gal et ail (1201^ develop a new 
model that could explain the formation of these hydrides 
in the gas phase while some of the modelers mentioned 
above claim surface formation. Concerning oxygen, one 
of the conclusions of deep ob servations of O2 towards 
NGC 1333 IRAS 4A protostar (|Yildiz et al.ll2013li is that 


the observed low molecular oxygen abundance is due to 
the freeze-out of atomic O onto grains. All these obser¬ 
vations lead to the conclusion that on the grains in dense 
clouds surrounding protostellar cores, there can be NH3 
and O present. 

On the experimental side, two important facts are of 
interest. First, conc erning NH3, it has been experimen¬ 
tally demonstrated (|Hidaka et al.l[2011r) that the succes¬ 
sive hydrogenation of N atoms trapped in a solid N2 
matrix at 10 K leads to NH3 formation. Second, the 
formation of NH2OH has been observed after irradiat- 
ing an ammonia-w ater ice at 90-130 K with UV photons 
(iNishi et al.lll98 4fl. and at 10 and 50 K with 5 keV elec¬ 
trons ( Zheng fc Raised l2010fl . The chemi cal processing: 
of H2 O-NH3 ices induced by heavy ions (|Bordalo et al.l 
[Ml leads also, by radiolysis, to the formation of 


N2O, NO, NO2 and NH2OH. As mentioned above, 
(jBlagoievic et al.ll200l showed that the reaction of hy¬ 


droxy lamine ions, NH2OH+, with acetic or propanoic 
acid makes ionized glycine and alanine. Much earlier, 
from about 1960 to mid 1980s, the NH3-I-O reaction has 
been the sub ject of numerous experiments reviewed by 
(jGohenlflQSTti . They have been performed at moderate 
temp erature (450 to 850 K~) in flow ing and static sys¬ 
tems (|Perrvlll9^ : lBaulch et al.lll984[l where ground state 
0(^P) is obtained by laser photolysis , or at higher tem¬ 
perat ures (850 to 2200 K) in flames (iFenimore fc ,lon^ 
Il961t) or in shock tube experiments (|Fuiii et al.lll98l 7 
Even if such a reaction has never been studied at very low 
temperature, it has been proposed that excited NH3O in- 
termedia te could rearrange t o give stable hydroxylamine 
NH2OH (|Baulch et al.lll98l . The formation of complex 
organic molecules involving atomic addition reactions on 
the cold surface of ices-covered interstellar grains in di¬ 
verse e nvironments, has been exten sively reviewed re¬ 
cently ([Herbst fc van Dishoedd [200911 . As NO has been 
detected in the gas phase towards a few high-mass star¬ 
forming regions, it has been suggested that hydroxy¬ 
lamine could be formed in the gas phase through non- 
energeti c hydrogenation of NO ice under dark cloud con¬ 
ditions (jCharnlev et al.ir2001f) . But recently, an NO ice 
hydrogenation surface reaction at low temperature has 
shown that an efficient route to NH2OH f ormation was 
also po ss ible with interstellar r elevant ices ([Gongiu et al.l 
l2012alibl : iFedoseev et al.ll?012h . It takes place at 10 to 
15 K at both NO submonolayer and multilayer cover¬ 


age. Such a process occurs without external energetic 
input but implies several successive steps of hydrogena¬ 
tion. It is worth mentioning a crossed molec ular beam 
study of the 0(^D)-|-NH3 reaction at 295 K (jShu et al.l 
I2OOII) forming vibrationally excited NH2OH with a sub¬ 
sequent dissociation into two different reaction channels: 
OH-fNH2 and NHOH/NHaO-fH with 90% and 10% yield 
respectiv ely. This reaction has been also studied theo¬ 
retically (iWang et al.ll2004ll . Finally, very recently, hy¬ 
droxylamine have been formed in a two-step mechanism 
through the reaction of a mmonia with hydro xyl radicals 
in a Neon matrix at 4 K (|Zins fc Krimll20T^ . 

However, most of the objects mentioned above, if not 
all, either in high-mass or low-mass star-forming regions 
as well as in proto-stellar cores, share the common char¬ 
acteristic of having a large abundance of NH3 (and also 
H2O), in any case far higher than the NO abundance. 
This latter molecule appears to be present in a restricted 
number of objects relating to high-mass star-forming re¬ 
gions. The abundance of NH3 is a sserted in observations 
of low-mass young stellar objects ([Bottinelli et al.ll2010H 
and even dark clouds (a s the Bok globule Barnard 68 
(|Di Francesco et ani2002tP and comets (as very recently 
detected on 67P/C-G by the ROSINA instrument aboard 
ROSETTA (|Altwegdl2014f) b Then the formation pro¬ 
cess in star-forming regions of hydroxylamine we present 
in this paper is a very likely more direct (single step) 
and realistic (considering the species inv olved) mecha¬ 
nism, compared to the one already known (jGongiii et al.l 
I2m2bll . 

Actually this process is well suited with respect to to 
hot-cores and protostellar cores. In addition to N result¬ 
ing from N2 dissociation under cosmic-rays irradiation, 
atomic O is also present for the same reason. Both N and 
O stick on the colder bare or ices covered grains, the ones 
deep into the cloud surrounding the core (A„ ~ 15 - 20). 
A very simplistic scheme could then be the following: (i) 
Very cold grains are covered by atoms (H, O, and N), ei¬ 
ther on the surface or in bulk, (ii) when the temperature 
begins to increase (>5-10 K), H atoms become mobile 
and form NH3 (and H2O) by hydrogenation surface re¬ 
actions , (iii) as the temperature continues to increase 
(>14 K in interstellar space conditions - corresponding 
in laboratory experiments to >20-25 K)), O atoms be¬ 
come mobile to finally react with NH3 (and H2O) already 
formed to to produce in particular NII2OH through an 
insertion mechanism. 

In this paper we present experimental evidence of the 
formation process of hydroxylamine via oxidation of NH3 
on dust grain analogs (Section 3). The experimental 
findings are then used in a 3-phase gas-grain model of 
a dark cloud (Section 4). A comparison of NH2OH for¬ 
mation mechanisms and a discussion about suitable envi¬ 
ronments where NH2OH could be observed are presented 
in Sections 5 and 6. 

2. EXPERIMENTAL SETUP 
2.1. Apparatus 

The apparatus was described elsewhere (|He 
iJing et ^l2013ll : here we summarize the main features 
that are important for this study. The experiments 
were carried out in an ultra-high vacuum main cham¬ 
ber connected to an atomic/molecular beam line. The 
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main chamber is pumped by a combination of a cryop- 
ump, turbomolecular pump, and ion pump and can reach 
1.2 X 10“^° torr after a bake-out. At the center of the 
chamber there is a 1 ^m thick amorphous silicate thin 
film sample grown on a 1 cm^ gold plated copper disk by 
electron beam physical vapor deposition. The detailed 
preparati on and character ization of the sample can be 
found in iJing et ahl (|2013ll . The sample can be cooled 
to 8 K by a liquid helium-cooled sample holder and be 
heated up to 450 K by a cartridge heater behind the sam¬ 
ple. At the beginning of each day, the sample is heated 
up to 400 K to clean it before cooling it down. The cham¬ 
ber pressure after cooling down is about 5 x 10“^^ torr. 
During exposure of the sample to the atomic/molecular 
beam, the pressure increases to about 1 x 10 “^° torr. 
At this pressure, the amount of water that stick on the 
sample from the chamber background gas is negligible. 
A triple-pass quadrupole mass spectrometer (QMS) is 
mounted on a rotary platform to record desorbed species 
from the sample surface or to measure the beam com¬ 
position. The QMS is placed in a differentially pumped 
enclosure and is fitted with a cone aimed at a sample. 
The distance between the tip of the cone and the sample 
can be changed by repositioning the sample via a XYZ 
sample holder manipulator. This arrangement prevents 
products loosely adsorbed on sample holder parts from 
reaching the QMS, thus reducing the background level 
and improving the signal/noise. A radio-frequency dis¬ 
sociation source is mounted in the first stage of the triple 
differentially pumped beam line. A mass flow controller 
(Alicat MCS-5SCCM) is used to ensure stable gas flow. 

2 . 2 . Beam flux calibration 

The beam fluxes are calibrated by TPD experiments. 
With the the amorphous silicate sample kept at 70 K, 
ammonia gas is introduced from the molecular beam line. 
The gas flow is set to 0.3 seem on the mass flow con¬ 
troller. After ammonia exposure, the sample is cooled 
down to below 40 K and then heated up at 1 K/s to 
about 200 K to do a TPD. In both ammonia exposure 
stage and TPD stage, the gas phase molecules are mea¬ 
sured using the QMS. Since ammonia fragments in the 
QMS ionizer, both mass 17 amu (NH;j') and 16 amu 
(NH^) are measured. The fragmentation of ammonia 
into NH+/N+ (mass 15/14 amu) is found to be negli¬ 
gible by direct beam measurements. A series of TPD 
experimental runs with different ammonia doses are per¬ 
formed. The TPD spectra are shown in Figure [1] We 
calibrate the NH 3 coverage by analyzing the shape of 
the TPD profile. In the submonolayer region the TPD 
peak temperature should shift to lower values, because 
in the low coverage region molecules occupy the deep 
sites (higher desorption energies). Here a layer is actu¬ 
ally an equivalent layer, since it is unclear whether am¬ 
monia would form clusters or islands on the surface; little 
experimental evidence support the existence of clusters 
or islands. As the deep sites get filled up with increasing 
exposure, the TPD pr ofile and peak tem perature shift 
to lower temperatures (jHe fc Vidalill2014[l . In the mul¬ 
tilayer region , the TPD spec tra should show a common 
leading edge (|Kolasinkill2008ll . which is a typical first or¬ 
der desorption behavior. The trend of TPD traces shown 
in Figure [T] is in agreement with the one obtained by 



Figure 1. TPDs with different exposure of NH 3 at 70 K. The 
exposure times are, from bottom to top, 0.5, 1, 2, 4, and 8 min¬ 
utes, respectively. Desorption energy distribution of NH 3 calcu¬ 
lated from the 2 minutes exposure curve (1 ML deposition) is in 
the inset. 

iBolina et al.l (|2005ll for NH 3 desorption from a graphite 
surface. With the current experimental settings, 1 ML 
NH 3 coverage is achieved by 2 minutes exposure with 
0.3 seem gas flow. Thus the NH 3 flow is 0.5 ML/minute, 
or equivalently, 5 x 10^^ cm“^minute“^, assuming 1 ML 
^ 10^^ cm~^ . Foll owing the procedures as described in 
iHe fc Vidalil (j2014ll . the desorption energy distribution 
of NH 3 from amorphous silicate surface can be obtained 
by direct inversion of the 1 ML trace in Figure [T] The 
resulted distribution is shown in the inset of Figure [TJ 
Additional TPD runs are performed at a surface tem¬ 
perature of 10 K, 30 K, and 50 K. They show almost 
identical TPD curves as the one at 70 K. This suggests 
that the sticking of ammonia on the silicate surface at 70 
K is unity, and the desorption rate at 70 K is negligible. 

In the atomic oxygen exposure, the O 2 gas flow pass¬ 
ing through the flow controller is set to 0.1 seem. At 
this flow the dissociation rate of O 2 is measured to be 
42%. The calibration of oxygen flux needs to be done 
differently from ammonia because O 2 is volatile and 
the sticking is not necessarily unity. The direct beam 
intensities of O 2 gas at both 0.3 seem and 0.1 seem 
are measured. Since it is known that 0.3 seem corre¬ 
sponds to 0.5 ML/minute, the flow at 0.1 seem can be 
calculated correspondingly. The O 2 flow at 0.1 seem 
is about 0.34 ML/minute (3.4 x 10^^ cm“^minute“^). 
When the RF is on, the beam intensities of O and O 2 
are 0.29 ML/minute (2.9 x 10^"^ cm“^minute“^) and 
0.20 ML/minute (2.0x10^^ cm“^minute“^), respectively. 
With the radio frequency (RF) power on and the oxygen 
sent into the dissociation source, beam contamination is 
checked. The main contaminant is NO (mass 30 amu) 
due the small leak of air into the dissociation source. 
The NO (mass 30) signal is less than 3% of the O sig¬ 
nal, which is trivial in the context of the experiments 
performed here. 

2.3. Experimental procedures 

The ammonia oxidation reactions were studied in am¬ 
monia and atomic oxygen sequential exposure TPD ex¬ 
periments. The surface was covered with ammonia be¬ 
fore introducing atomic oxygen. Two ammonia coverages 
were used, 2 ML and 1/4 ML, representative of multilayer 
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Figure 2. A typical TPD of a NH 3 + O sequential exposure 
experiment at 70 K in which the QMS records simultaneously mul¬ 
tiple signals. The temperature ramp (black line) shows that after 
the exposure is terminated, the sample is cooled to 40 K before the 
heating is started. 

coverages and submonolayer coverages. After prepara¬ 
tion of the ammonia sample, the residual ammonia in 
the beam line and dissociation source was cleaned by 
flushing the beam line several times using oxygen. This 
ensures that almost no ammonia was mixed with oxy¬ 
gen. The sample temperature in the exposure stage was 
chosen to be 70 K so that O/O 2 /O 3 does not stick onto 
the surface while the sticking of NH 3 is unity. After the 
exposure stage, the surface was cooled down to below 40 
K and then heated up linearly at 1 K/s to above 320 K 
to desorb species from the surface. The QMS recorded 
simultaneously the signals of various masses during both 
exposure and TPD stages. Figure [2] shows the TPD spec¬ 
tra of a typical experimental run. 

3. RESULTS AND ANALYSIS 
3.1. Multilayer NH 3 + O 

In the NH 3 + O sequential exposure TPD experiments, 
mass 16, 17, 18, 30, 32, and 33 amu were measured by 
the QMS, see Figure [2l Mass 16 amu is due to the frag¬ 
mentation of NH 3 in the QMS ionizer; mass 17 amu is 
due to NH 3 and fragmentation of H 2 O; mass 32 amu is 
due to O 2 and fragments of O 3 . Mass 30 amu could be 
due to NO from the beam line (but the amount is small), 
or fragments of NH 2 OH; mass 33 amu could be due to 
NH 2 OH. HO 2 has the same mass but it is unlikely to 
be there because of the lack of detection of the prod¬ 
ucts it fragments into. The mass 33 amu peak centered 
at around 180 K is accompanied by mass 30 amu and 
very small signals with masses 16, 17 and 32 amu at the 
same temperature; this suggests that the mass 33 amu 
peak at 180 K is due to NH 2 OH. This agrees with the 
TPD peak attribution in ICongiu et al.l (l2012aD ( in the 
range of 160-200 K, peaking at T~ 190 K) obtained in 
an NO hydrogenation experiment and confirmed by in¬ 
frared measurements (RAIRS). Mass 18 amu is due to 
water. 

Figure [3] shows the mass 17 amu (NII 3 ) signal after 
depositing various doses of O on top of 2 ML (4 min¬ 
utes exposure) NH 3 . As the O exposure time increases, 
the NH 3 peak decreases and the peak position also shifts 
to higher temperatures. This is because the top layer of 
NH 3 is gradually converted to NH 2 OH or other products, 
thus hindering the desorption of NH 3 molecules under¬ 



Figure 3. Mass 17 amu QMS signal of the TPD from an ice 
consisting of NH 3 (4 min or 2 ML) followed by exposure to various 
doses of O at 70 K. 

neath. After about 6 minutes of O exposure, the first 
peak almost disappears, indicating the top layer of NH 3 
is almost all converted to NH 2 OH or other products. 

In Figure[4]two desorption peaks, peak A and peak B, 
are visible. Peak A shows up at O exposure as low as 
0.5 minutes. This indicates that the ammonia oxidation 
reaction is efficient. As the O exposure increases from 0 
to 4 minutes, the area of peak A increases at first, but 
then it decreases from 4 minutes to 16 minutes of O ex¬ 
posure. Peak A is followed, at larger O exposures, by 
peak B at 260 K. Notice that peak B starts to ap¬ 
pear at 4 minutes of O exposure time at the expense of 
peak A. As mentioned above, peak A falls in the temper- 
ature range of the T PD mass 33 amu peak obtained by 
ICongiu et al.l ()2ni2aD . Therefore we attribute peak A to 
the desorption of NH 2 OH. Concerning peak B, however, 
in a s ubsequent paper on an NO 2 hydrogenation experi¬ 
ment (lIoDDolo et al.ll20l1l , a mass 33 peak, confirmed by 
infrared measurement spectra, appears at a temperature 
T^ 245 K, close to the temperature of our second peak 
(peak B in our Figure |4|). Nevertheless, oddly there is no 
mention in their paper of a peak A at 160 K. 

To find out what peak B represents, we show in Figure 
[5] a comparison of TPD traces for various masses for 2 
ML of NH 3 -b 8 or 16 minutes of O exposure, see left and 
right panels of Figure [5l respectively. Peak A and peak 
B are marked by vertical lines. In both panels of Figure 
O peak A (see trace of mass 33 amu) is accompanied by 
desorption of mass 30 amu, 17 amu, and 16 amu, while 
this is not true for peak B. This suggests that peak B is 
not due to NH 2 OH. I t could be the product of a frag¬ 
mentation of a dimer (|Del Benelll972h or of an oxidation 
product of NH 2 OH. At about 280 K, there is a peak C 
showing up for mass 30 amu, 17 amu, and 16 amu, but 
only when the O dose is high (16 minutes of O expo¬ 
sure), see right panel of Figure 111 It could be due to a 
yet unidentified product formed in a further oxidation of 
NH 2 OH. Peak B and peak C differs in both position and 
shape, therefore they should be attributed to different 
species. 

3.2. Submonolayer NH 3 + O 

Experiments were also carried out at submonolayer 
NH 3 coverage. We exposed the silicate sample to NH 3 
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for 0.5 minutes, or about 1/4 of a layer, with the same 
procedures as the ones followed to obtain the data in Fig- 
ure[Tl This was followed by exposure to O. Mass 16 amu 
was chosen to represent the NH 3 amount because the 
mass 17 amu signal has a non-trivial contribution from 
water fragmentation when the NH 3 signal is weak. In 
Figure[ 6 ]we show the integrated mass 16 amu TPD trace 
for different exposures to O. With O exposure from 0 to 
2 minutes, the amount of NH 3 follows more or less an 
exponential decay. This is because when the O amount 
is small, oxidation dominates and the NH 3 destruction 
rate is proportional to the NH 3 amount. With a further 
increase in O exposure, the NH 3 decay does not follow 
a simple exponential decay anymore because of possible 
secondary reactions. In Figure [ 6 l a straight line is fit¬ 
ted to the loge plot from 0 to 2 minutes. The slope is 
—0.46 ± 0.06 minute”^. In an exponential decay, the 
amount of NH 3 on the surface should follow exp(—uf/t), 
where cr and cj) are the reaction cross-section area and 
the beam flux, respectively. From it, we obtain the cross 
section a = —slope/(f> = ( 1.6 ± 0 . 2 ) x 10 “^^ cm^. 

3.3. Reaction energy barrier of NH^ + O 

In typical ISM conditions, it is unlikely that situations 
arise similar to the ones that led to the appearance of the 
second TPD peak of mass 33 amu at high temperature 
attributed to a subsequent oxidation process. Therefore 
we focus on the first oxidation: 

NH 3 + O ^ NH 2 OH (I) 

When the surface is fully covered with NH 3 and the 
surface temperature is 70 K, O could react with NH 3 
also via the Eley-Rideal or hot-atom mechanisms, i.e. 
reaction without complete thermal accommo dation of 
the o xygen atom. Based on prior experiments (jHe et al.l 
[ 2 ?fn .‘ we assume that at 70 K the residence time of O 
on NH 3 ice is negligible. We also assume that all the 
reactions take place during the exposure stage instead of 
the TPD stage, because at 70 K the mobility of oxygen 
should be relatively high. The degree of conversion of 
NH 3 can be gauged by looking at what remains of the 
NH 3 layers after various doses of oxygen. Figure [3] shows 
the TPD mass 17 signal (NH 3 ) for various doses of O. 
The NH 3 dose is fixed at 4 minutes, which is equivalent 
to about 2 ML. As the O exposure time increases, the 
NH 3 peak decreases, and the peak position also shifts 
to higher temperatures. This is because the top layer of 


NH 3 is gradually converted to NH 2 OH or other products, 
thus preventing NH 3 molecules underneath from desorb¬ 
ing. After about 6 minutes of O exposure, the first peak 
almost disappears, indicating that the top layer of NH 3 is 
almost all converted to NH 2 OH or other products. Thus 
it takes at most 6 min of O (1.74 ML) to convert 1 ML 
of NH 3 , suggestion an efficiency for the reaction of each 
oxygen atom of at least 0.575 (=1/1.74). 

Using this efficiency, an upper limit for the activation 
energy barrier, Ea, for the O -I- NH 3 reaction may be 
estimated. 

We firstly assume that the incoming oxygen atom is 
thermalized on the surface, before overcoming the acti¬ 
vation energy barrier to reaction. This thermalization 
is the most likely immediate outcome of the impact, ex¬ 
cept in a small minority of ’’direct-hit” trajectories, as 
the oxygen atom dissipates energy as it is drawn into the 
multi-dimensional potential of the surface. The atom is 
nevertheless available for immediate reaction if the ac¬ 
tivation energy barrier can then be overcome; in such a 
case, the reaction may be classed as Eley-Rideal, as it is 
not mediated by a diffusion process. 

Assuming that the main loss mechanism for an oxy¬ 
gen atom on a pure NH 3 surface is either reaction with 
NH 3 or thermal desorption, the reaction efficiency may 
be formulated thus: 

1/1.74 = exp[-F;A/r]/(exp[-AA/T]+exp[-Ades(0)/r]) 

( 2 ) 

where E^esiO) is the desorption energy for an oxy¬ 
gen atom and T is the surface temperature. Adopting 
Edes(O) = 1500 K and T = 70 K, this expression pro¬ 
duces a value i?A( 0 -l-NH 3 ) = 1479 K. 

This value constitutes an upper boundary of the reac¬ 
tion energy barrier because we assume a sticking coeffi¬ 
cient of unity, and thence that each oxygen atom that 
reaches the surface participates in a reaction. If this is 
not the case, then the reaction energy barrier for the first 
oxidation should be even smaller. 

The above treatment does not consider the possibil¬ 
ity that surface oxygen atoms could meet each other 
and react, reducing the efficiency of the reaction. Un¬ 
der conditions in which such reactions could be impor¬ 
tant, then in order for the O -I- NH 3 reaction to be ef¬ 
ficient and the 0 - 1-0 reaction to be minimized, the 
reaction barrier for O -I- NH 3 would need to be lower 
than the oxygen diffusion barrier. This quantity has not 
been measured in the laboratory, but there are rules- 
of-thumb based on experimental data from other sys¬ 
tems. The ratio of the diffusion energy barrier to the 
desorption ener gy is typically foun d to be '^0.3 for crys¬ 
talline surfaces l)Bruch et al.ll2007ll . For rou gh surfaces, 
values between 0.5 and 0.8 h ave been used (IKatz et al. 
II 999 I: ^rrod &: Herbstll200^ . although iGarrod fc Paulv 
(j2011l l found that values lower than ~0.4 were most con¬ 
sistent with observed interst ellar ext i nction thresholds 
for H 2 O, CO and CO 2 ices. iGarrodI (j2013[ l adopted a 
value of 0.35, which would yield an upper limit on the O 
-|- NH 3 activation energy barrier as low as 525 K. 

3.4. Control experiments 

Control experiments were carried out to verify whether 
ammonia react with molecular oxygen or ozone. We ex- 
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O exposure time (min) 

Figure 6 . Integrated TPD signs mass 16 amu from a deposit of 
NH 3 (0.5 minutes or 1/4 of a ML) + various exposures of O at 70 
K. 

posed the ammonia layers to O 2 and found no mass 33 
or 30 amu peaks in the TPD, signifying that there has 
been no conversion of ammonia to hydroxylamine in the 
presence of molecular oxygen. We also checked the re¬ 
activity of ozone with ammonia. Ozone was prepared 
on a clean silic ate following the procedure described in 
I He et al] l)2014r) ; then, 2 ML of ammonia were deposited 
on it. Again, there was no mass 33 or 30 amu peaks in 
the following TPD, see Figure El 

4. MODELING 

In order to test the astrophysical importance of the 
O -f NH 3 —>■ NH 2 OH reaction investigated here, we in¬ 
corporate this surface mechanism into the recent three- 
phase gas-grain chemical kinetics model MAGICKAL 
(Garrod 2013). The model parameters used here are in¬ 
tended to approximate those present under dark cloud 
conditions, during which significant ice mantles are 
formed on the surfaces of dust grains. The model uti¬ 
lizes the chem ical network a nd initial elemental abun¬ 
dances used bv iGarrodI (j2013[ l. assuming a generic, static 
dark-cloud gas density of nn = 2 x 10 ^ cm“^ and vi¬ 
sual extinction Ay = 10. A gas temperature of 10 K 
is used in all model runs. We assume a binding energy 
for atomic oxygen of 1500 K, based on an experiment 
that measured the des orption energy of O from compact 
amorphous water ice (iHe fc Vidaliil2014ll . Diffusion bar¬ 
riers are set to a uniform fraction of binding energies for 




Figure 7. Mass 33 QMS signal in the TPD from NH 3 -|- O 2 and 
O 3 + NH 3 deposited at 70K compared with that of NH 3 + O at 
50 K. Curves are offset for clarity. 

all species, EdiS = 0.35 Edes, following (jGarrodI 1201311 . 
Due to the low temperatures at which grain-surface ice 
mantles are formed, only surface chemistry is switched on 
in this model, although the three-phase treatment con¬ 
siders separate grain/ice surface and bulk mantle popu¬ 
lations for all species (as well as gas-phase abundances). 
The three-phase model allows the composition of each 
layer within the ice mantle to be traced as it is de- 
posited during t h e chem ical evolution of the cloud (see 
iGarrod &: Paulvl (|2011ll l. Only Langmuir-Hinshelwood 
(i.e. diffusional) processes are considered in this model, 
as the Eley-Rideal mechanism tends to require surface 
coverages of close to unity (in models of astrophysical 
grains) to compete with L-H mediated reactions. 

Reactions are also included in the network allowing 
NO molecules to be hydrogenated to NH 2 OH by the se - 
quential addition of atomic H. Following (lGarrodll2013[l . 
this includes a barrier-mediated reaction, H -|- HNO —>■ 
HNOH, as well as an alternative H-abstraction branch, 
H -b HNO —>■ H 2 -b NO, for which the b arrier is assumed 
to be 1500 K in both cases. Following (jHaseeawa et al.l 
Il992f) . tunneling through activation energy barriers is 
treated using a rectangular barrier treatment, with a de¬ 
fault width of 1 A. The faster of tunneling and thermal 
penetration is used, according to the grain temperature. 

The model described above is used to produce a grid 
of models for a range of dust-grain temperatures from 10 
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- 20 K, and using activation energy barriers, i?Aj for the 
O + NH 3 —NH 2 OH reaction ranging from 0 - 2500 K. 
A smaller set of models is also run to compare the effects 
of NH 2 OH produced by this mechanism, with cases in 
which the barrier to the H + HNO —?► HNOH reaction is 
altered. 

4.1. Model results 

FigurelSl panel (a), shows the composition of each layer 
in the grain mantle ice for the model run corresponding 
to i?A=0 and Tg^ain = 10 K. Typically observed com¬ 
ponents are shown, along with NH 2 OH. In this case, 
NH 2 OH production is relatively small, and all such pro¬ 
duction is a result of the hydrogenation of NO. At 10 
K, oxygen atoms are insufficiently mobile to be able to 
reach and react with NH 3 before they are hydrogenated 
to produce water, even under the assumption of a zero 
barrier for the O -I- NH 3 reaction. Here, the majority 
of NH 2 OH is produced via NO hydrogenation, which is 
strongly peaked at late time/outer layers. 

Panels (b) and (c) show the results at 14 K, for acti¬ 
vation energy barriers of Aa=1000 and 2000 K, respec¬ 
tively. At this temperature, oxygen atoms are sufficiently 
mobile to meet ammonia molecules on the surface be¬ 
fore being hydrogenated. In the i?A=1000 K case, the 
reaction barrier is low enough (in comparison to the dif¬ 
fusion barrier for the oxygen atom to move away from 
the ammonia molecule) to allow reaction to occur with 
high efficiency. In this case, the majority of NH 2 OH is 
formed through this reaction, and NH 2 OH abundance 
closely follows that of NH 3 , never falling below a factor 
of ten lower than NH 3 at any layer in the ice. NH 2 OH 
maintains a fairly steady abundance throughout the ice, 
up to around 1 Myr (170 ML), but the NH 2 OH produc¬ 
tion is enhanced at all times/depths. In the Aa=2000 K 
case, the reaction barrier is too high to allow reaction be¬ 
fore hydrogenation of oxygen occurs, making the NH 3 -|- 
O reaction inefficient; here NH 2 OH production is again 
the result of NO hydrogenation. 

Figure [9] shows the final (maximum) grain-mantle 
abundances of NH 2 OH as a function of total hydrogen 
in the cloud, for model runs assuming a range of grain 
temperatures from 10 - 20 K, and activation energies for 
the O -I- NH 3 reaction ranging from 0 - 1500 K. Results 
for the Ea= 500 K case are not shown, as they are identi¬ 
cal to the Ea = 0 case; both values are sufficiently below 
the atomic oxygen diffusion barrier assumed in the model 
(525 K), such that the reaction with NH 3 occurs at an ef¬ 
ficiency close to unity, if the temperature is high enough 
to allow the oxygen atom to diffuse and reach an NH 3 
before it is hydrogenated and becomes immobile. This 
is the case for models of temperature 14 K and above; 
all models at 12 K or less produce the same results, with 
NH 2 OH produced by NO hydrogenation. 

Assuming a barrier of 1000 K, the reaction is still rel¬ 
atively efficient, falling by a factor ~2 at grain tempera¬ 
tures of ^^16 K and above, compared to the Ea =0/500 
K results. However, with a 1500 K barrier, reaction is 
inefficient, and yet higher barriers produce the same re¬ 
sults; in such cases, NH 2 OH is produced overwhelmingly 
by NO hydrogenation. A run with the O -I- NH 3 bar¬ 
rier set to 1250 K allows this reaction to become slightly 
dominant over the NO hydrogenation mechanism, pro¬ 
ducing a small, though significant, increase in NH 2 OH 


production. 

Models were also run to compare the case where the 
hydrogenation of NO by H atoms is assumed to occur ei¬ 
ther (i) without barriers or (ii) not to occur at all. Case 
(i), removal of the barrier to H -|- HNO —>■ HNOH (while 
the alternative hydrogen-abstraction branch retains its 
barrier of 1500 K), in fact results in no appreciable dif¬ 
ference from the runs that assume a barrier of 1500 K for 
H -I- HNO (the zero-barrier results are thus not plotted 
explicitly in the figure). NH 2 OH production through this 
mechanism is therefore already at its highest efficiency. 

The alternative regime, case (ii), in which the H -|- 
HNO —>• HNOH barrier is increased to arbitrary height, 
such that the reaction does not occur, results in negligi¬ 
ble NH 2 OH production at 12 K and below, because all 
NH 2 OH formation is dependent on the O -t- NH 3 reac¬ 
tion, which is inefficient at these temperatures. Above 12 
K, NH 2 OH is formed efficiently via O -|- NH 3 , and there 
is only a small difference between these results and those 
in which the H -|- HNO reaction is allowed to occur at 
maximum efficiency. 

5. DISCUSSION 

Experiments of oxygen exposure on submonolayer and 
multilayer ammonia ices show an efficient formation of 
NH 2 OH with a reaction energy barrier that can be as low 
as ^ 52577 and no higher than 1479K. The experiments 
were done on samples at 70 K to avoid contamination of 
molecular oxygen and ozone. At this temperature NH 3 
is stuck on the surface while O has a short residence 
time. Experiments have been performed at both bilayer 
and submonolayer coverage. Simulations were done to 
find what rates of hydroxylamine formation one would 
obtain in simulated dense cloud conditions. According 
to the chemical models, at dust temperatures around 14 
K or higher, the O -I- NH 3 reaction is found to be an 
efficient mechanism for the production of NH 2 OH. While 
temperatures this high are not the typically considered 
dark-cloud dust temperatures, it is certain that the dust 
in dark clouds must pass through such temperatures on 
their way to the 8 - 10 K that is more typically assumed. 
Furthermore, interstellar dust will again pass through the 
^14 K threshold at which this reaction becomes efficient 
as part of the evolution of hot cores. 

6 . ASTROPHYSICAL IMPLICATIONS 

NH 3 and H 2 O have been detected in the gas phase 
around transient protostellar cores or inferred for O from 
the strong O 2 depletion. These cores are also classi¬ 
fied as hot cores or hot corinos in the case of mas¬ 
sive and medium/low mass nascent stars, respectively. 
These cores whose temperature can be as high as 100 
K close to the protostellar source (and even 250 K for 
the densest ones) are much colder far from the core due 
to the strong extinction provided by the dusty molec¬ 
ular clumps in which they are embedded. A possible 
scenario for the coldest regions where NH 3 and O are 
stuck on the grains and approaching the core is that in 
a first stage (at temperature in the 14-50 K range) O is 
released in the gas phase and is then able to collide with 
NH 3 covered grains, inducing chemical reactions that can 
synthesize new species. These would later desorb from 
the grains when the core temperature increases. Alterna¬ 
tively, NH 2 OH may be produced at lower temperatures. 
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Figure 8 . Fractional composition of ices within each layer as a 
function of depth into the ice, lower horizontal scale, or, equiva¬ 
lently, as a function of age of the cloud, top horizontal scale. Panel 
(a) The activation energy barrier for the O-I-NH 3 reaction is 0 K 
and the temperature of the ice is 10 K. Panel (b) F'a= 1000 K and 
T=14 K. Panel (c) F;a= 2000 K and T=14 K. 



Figure 9. Final model abundances of ice-mantle NH 2 OH with 
respect to total hydrogen. The grain temperature is varied for a 
selection of activation-energy values for the NH 3 -|- O reaction. 
The dashed line shows the results assuming no barrier for either 
the NH 3 + ^ NH 2 OH or H -h HNOH ^ HNOH reactions. 


during the formation of the ice mantles, to be released 
at later times. (Models of such scenarios are beyond the 
scope of this primarily experimental study, but will be 
considered in future work). 

Consequently, the present laboratory astrophysics ex¬ 
periment on O interaction with NH 3 should be useful to 
interpret the observational data. Then the question is of 
the non-observation of hydroxylamine while, at the same 
time, it has been demonstrated that it can be made in 
the laboratory in conditions almost similar to those in 
protostellar cores.The upper limit relative abundance of 
NH^ OH with respect to H^. as deduced from observa¬ 
tions (iPulliam et al.l[2C)12ll is found to be in the 3 x 10“® 
to 8 X 10“^^ range, depending on the object. Interest¬ 
ingly, combining this result with the relative abundance 
of NH 3 with respec t to H 2 in a high-mas s star-forming 
region, 4.6 x 10“® (iBattersbv et al.l[201411 (a value of 2 
X 10“® is already reported in th e case of a cold dense 
cloud like TMC-1 (I Turned [200011 '). we deduce an aver¬ 
age fractional abundance of NH 2 OH with respect to NH 3 
in the 0.7 x 10“^ to 2 x 10“^ range, close to the val¬ 
ues deduced from the proposed model. Several reasons 
may explain the non-observation: the yield of such re¬ 
actions could be quite low due to the narrow tempera¬ 
ture range in which they have to take place as well as 
a possible spatial confinement of dusty clouds where the 
required temperature conditions for f ormation are ful¬ 
filled. In addition, recart ob servations (jSakai et al.ll2013l 
l2014airhl : iSanna et al.ll2014ll have shown that the orien¬ 
tation of the protostellar core in the plane of the sky 
could be also important, particularly if it is associated 
to infall-disk-outflow systems and rotation. It could be 
expected that the amount of chemical species resulting 
from the NH 3 -I-O interaction will be almost equal both 
in high-mass and medium- to low-mass stars because of 
a higher species density in the first case and a longer in¬ 
teraction time in the other. All these conditions could 
explain the difficulty to observe such species as hydrox¬ 
ylamine that have, however, been made in the labora¬ 
tory. ALMA could provide the required sensitivity and 
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the sub-arcsecond a ngular resolution to observe them. 

A recent paper (|Cleeves et al.l l2014f) demonstrates, 
based on a physical and chemical model of the proto¬ 
planetary disk of the solar system, that the relative 
deuterated water abundance observed on Earth cannot 
be explained by its formation in the disk itself but should 
results in part from a primeval enrichment by interstellar 
ices. All species embedded in these ices may then have 
survived and have been injected into the disk and later 
on in the planetary system. This discovery confirms the 
paramount importance of observing prebiotic species like 
glycine or even smaller building blocks of such molecules 
like hydroxylamine, which would emphasize the possible 
extraterrestrial origin of life. 
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